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We describe a new mechanism of tunneling between period-two vibrational states of a weakly non- 
linear parametrically modulated oscillator. The tunneling results from resonant transitions induced 
by the fast oscillating terms conventionally disregarded in the rotating wave approximation (RWA). 
The tunneling amplitude displays resonant peaks as a function of the modulation frequency; near 
the maxima it is exponentially larger than the RWA tunneling amplitude. 
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I. INTRODUCTION 

Many systems of current interest can be modeled by 
modulated nonlinear quantum oscillators. Examples 
range from Josephson junction based systems jT| to opti- 
cal cavity modes [2] , electrons in a Penning trap [5] , and 
opto- and nano- mechanical systems [4]. The oscillator 
dynamics is often characterized by well-separated time 
scales: the reciprocal eigenfrequency u!q and a much 
longer time related to the vibration decay and nonlinear- 
ity. A standard approach to the analysis of the dynam- 
ics is based on the rotating wave approximation (RWA), 
where one separates slow variables, like the vibration am- 
plitude and the slow part of the phase, and disregards the 
effect of fast oscillating terms on their evolution. 

An important quantum effect in modulated systems 
is dynamical tunneling |5j. ft can be understood for a 
parametric oscillator, which is excited by modulation at 
frequency ujp close to 2ojq. Classically, a weakly nonlinear 
oscillator can have two states of vibrations at frequency 
wp/2, which have the same amplitudes and differ in phase 
by 7r [BJ. Quantum fluctuations cause tunneling between 
these states [7J [8]. Similar tunneling, which should be 
distinguished from dissipative switching [9, 10| , is known 
also for other types of vibration bistability [TT] . 

In this paper, we show that the tunneling rate of 
a parametrically modulated oscillator can be exponen- 
tially increased by processes caused by the fast oscillating 
terms oc exp(±mwj?i), (n — 1,2,.. .) disregarded in the 
RWA. This happens where the difference of the appropri- 
ate eigenvalues of the RWA Hamiltonian becomes close 
to nhuip. The level configuration is of A- type. The two 
lowest RWA levels are degenerate (disregarding tunnel- 
ing), with the wave functions localized near the period- 
two vibrational states, whereas the upper-level state is 
delocalized, see Fig. 1. The three states are resonantly 
mixed by the fast-oscillating non-RWA terms. The asso- 
ciated breakdown of the RWA is a purely quantum effect 
with no classical counterpart. 

The tunneling enhancement we consider is somewhat 
reminiscent of photon-assisted tunneling from a poten- 




neling 

uisi = F/2uip 



FIG. 1. (a) The dimensionless RWA Hamiltonian g(Q,P), 
Eq. pi, for (i = 0.95. The minima of g(Q,P) correspond to 
the parametrically excited vibrational states, in the presence 
of weak dissipation, (b) The cross-section g(Q, P — 0) with 
a sketch of RWA quasienergy levels; the arrows indicate reso- 
nant transitions due to the fast-oscillating corrections to the 
RWA. Also indicated are the typical energy scales. 



tial well, which is now broadly used in quantum infor- 
mation processing |12J. There photon absorption reso- 
nantly accelerates tunneling decay if the photon energy 
hup coincides with the intrawell level spacing, since the 
decay rate of the excited state largely exceeds that of the 
ground state. In contrast to systems displaying photon- 
assisted tunneling, a parametric oscillator is bistable due 
to the modulation, which forms the very barrier for tun- 
neling in phase space. This leads to a different physics 
and requires a different description. 

We study moderately strong resonant modulation 
where the nonlinear part of the oscillator vibration en- 
ergy remains small compared to the harmonic part. This 
makes the oscillator different from modulated strongly 
nonlinear systems where much attention has attracted 
chaos-assisted 13J and nonlinear resonance- assisted tun- 
neling [13] . 

For the rate of tunneling between the vibrational states 
to be small, the effective RWA tunneling barrier AC/ 
should largely exceed the RWA level spacing hu> s \; fre- 
quency o> s i characterizes the oscillator dynamics in the 
rotating frame, oj s \ <C up rj 2u>q. The ratio AU/Hujf can 
be arbitrary, it does not emerge in the RWA. We will be 



interested in the case where AJ7 ~ huip. In this case the 
effect of resonant admixture of the RWA states by the 
non-RWA interaction, is most pronounced, see Fig. 1. 

For ujp 3> Wgi) the states resonantly mixed by the non- 
RWA terms in the oscillator Hamiltonian overlap very 
weakly. We develop a method that allows us to calculate 
the relevant exponentially small matrix elements and to 
show that they are nevertheless sufficiently large to lead 
to exponential resonant enhancement of tunneling. 



II. RWA QUASIENERGY STATES AND THEIR 
RESONANT MIXING BEYOND THE RWA 

The Hamiltonian of a parametrically modulated oscil- 
lator with coordinate q, and momentum p reads 



P 



1 



H(t) = y + ~gf[<4 +Fcos(u F t)] + ^q 4 



(1) 



We assume that the modulation amplitude F and the 
nonlinearity are comparatively small, F 1 j(q 2 ) <C U!q, 
and the modulation frequency up is close to resonance, 
\uip — 2ljq\ <C uq; for concreteness we set F, 7 > 0. Of 
primary interest is the range of the modulation parame- 
ters F and up where, in the presence of weak damping, 
the oscillator has two almost sinusoidal stable classical 
vibrational states with frequency ujp/2 [BJ. 

To study the oscillator dynamics, we switch to the 
rotating frame using a standard transformation U(t) — 
exp[— iojpaJ at/2] and introduce dimensionless slow vari- 
ables Q and P, WqU = C[Pcosui F t/2-Qsmu F t/2] 
and WpU = -{ui F /2)C[Psmui F t/2 + Qcosui F t/2}. 
Here a and a* are the ladder operators and C = 
(2F/37) 1 / 2 . The Hamiltonian in the rotating frame 
H = WHU- ihWU reads 

H=(F 2 /6 1 )[g(Q,P) + h(Q,P,t)]. (2) 

The dimensionless operator 

\(Q 2 + P 2 ) 2 + l 2 (i-v)P 2 -li 

is independent of time (5). In contrast, the operator 



9 = - A (Q 2 + P 2 ) 2 + hl- v)P 2 -Ui + m)Q 2 (3) 



h 1 (Q,P)e 



-iujpt 



h 2 (Q,P)e 



-2iujpt 



H.c. 



is fast oscillating; hi t 2 are fourth order polynomials in 
Q, P, they do not contain small parameters and are 
given explicitly in the Supplemental Material. Functions 
g(Q, P) and h(Q, P) are symmetric with respect to inver- 
sion (Q,P) — > (-Q, —P) due to the periodicity of H(t). 
They depend on a single dimensionless parameter /j,, 

(j, = [(wf/2) - uj a ] /u B i, w 8 i = F/2u F . (4) 

We disregard corrections ~ ui s \/u>p. 

The commutation relation for the dimensionless coor- 
dinate Q and momentum P is 



where A is the dimensionless Planck constant. We assume 
that A < 1. Then quantum fluctuations are small on 
average. 

From Eq. (pi), the Schrodinger equation in dimension- 
less time r = tu) s \ is iXd^ = (g + h) 1 ^. Since h is 
periodic in time, this equation has Floquet solutions 
* £ (t + Th) = exp(-?erh,/A)* e (r). They define the di- 
mensionless quasienergies e [t^ — 2nu s \/ujp <C 1]. 

In the RWA the fast oscillating term h is disregarded. 
Then operator H becomes time-independent. The di- 
mensionless Hamiltonian g(Q,P), Eq. (pi, is shown in 
Fig. [T] It is not a sum of the kinetic and potential en- 
ergy. For \/j,\ < 1, g(Q,P) has two symmetrically lo- 
cated minima, g m i n = — (1 + /i) 2 /4, and a saddle point, 
9s = 0. In the presence of weak dissipation, the min- 
ima correspond to the period-2 vibrational states in the 
laboratory frame, which have equal amplitude and op- 
posite phase. The barrier height between the states is 
A[/ = (F 2 /6 7 )(.g 5 -3mi„). 

The eigenvalues g m of g give dimensionless quasiener- 
gies e m in the RWA. For A«l each well of g(Q, P) in 
Fig. [l] contains many levels, oc 1/A. Because the wells are 
symmetric, the intrawell states are degenerate in the ne- 
glect of tunneling. With account taken of tunneling, the 
eigenstates ip n (Q) of g are even or odd in Q. The dimen- 
sionless RWA tunnel splitting <5<? between the lowest-^ 
states was considered earlier [TJ, [5] . It is exponentially 
small, j log (5<7o I oc 1/A, and Sgo oscillates with /i/A [8]. 



The oscillating term h in the Hamiltonian (12j) mixes 
RWA-eigenstates. For remote states the mixing is ex- 
ponentially weak. However, it may become important 
where huip is close to the distance between the RWA lev- 
els, as it provides a new route for interwell transitions. 
Consider state ip n above the barrier top with dimension- 
al (r) 

less quasienergy g n and the two lowest states ip and ip Q 
in the left and right wells of g(Q, P) with quasienergy go 
in the neglect of tunneling, see Fig. [TVb) . The dimension- 
less detuning between the interlevel distance and huiF 
is A = X~ 1 (g n — g ) — (u)p/u> & \). If |A| <C 1, transi- 
ts r) 
tions i/iq — > ip n are resonant. The matrix elements 

(V'nl^ilV'o ) and (V'nl^ilV'o ) are ec l ua l f° r a symmet- 
ric ipn(Q) or have opposite signs for an antisymmetric 
VVi(Q); we denote their absolute value by h rcs . To first 
order in h, the amplitudes of resonantly coupled states 
ip n (Q),il>o ' r (Q) oscillate at dimensionless frequencies 



[(A 2 + 8A- 2 / lr 2 cs ) 1 /2 T | A | ]/2 . 



(6) 



[Q,P]=iX, X = 3jh/(Fu F ) 



(5) 



From Eq. (J6j), v± « \/2h rea /X for good resonance, 
X\A\/h TCS <C 1. In the dispersive regime, X\A\/h ICS 3> 1, 
interwell oscillations are characterized by frequency V- rj 
2h 2 es /X 2 A. As we show, in both cases v_ can be expo- 
nentially larger than the dimensionless RWA tunneling 
frequency 5go/\, which was disregarded in Eq. wm. 



III. MATRIX ELEMENTS FOR REMOTE 
STATES 

The relevant matrix elements of h can be found using 
the WKB approximation, in the spirit of Ref. [TSl In- 
terestingly, taking advantage of the conformal property 
of classical trajectories for the effective Hamiltonian g, 
one can find both the exponent and the prefactor in the 
matrix elements. For the term oc hi in h we write 

/•CO 

(ip n \hi\ij )=2Re dQh+(Q), 



h+(Q) = iP+(Q)hiMQ), (7) 

where ipo(Q) is one of the two tunnel-split lowest-g states 
[the symmetric or antisymmetric combination of ip Q (Q) 

and ipQ (Q) ]; VVi(Q) nas the same parity as ipo. Here we 
provide results for underbarrier states, g n < g$ = 0; for 
g n > the analysis is similar, see Supplemental Material. 
Function ipn(Q) IS an eigenfunction of operator g such 
that Re[i/>+((3)] = ip n (Q) and, in the classically accessible 
region of the semiaxis Q > 0, see Fig. 1, 

4>+(Q) « c n (d P g n y 1/2 exp (iX^S^Q) + *tt/4) . (8) 

Here, S n (Q) = J a , ,P(Q',g n )dQ' is the mechanical 
action counted off from the right turning point ap(g n ) 
and P(Q,g) is the classical momentum, 



P(Q,g)^^-l-Q 2 + 2^g~^. 



(9) 



In Eq. dsj) dpg n is dpg calculated for P = P(Q,g n ) and 
c n = [Tp (g n )/2]~ 1 / 2 , where t p '(g) is the period of clas- 
sical vibrations with quasienergy g. 

We shift the integration path in Eq. (jh to the upper 
half-plane, contour C in Fig. [2] (a). On this contour [TS] 



/i + ss [c n cohi(Q,-P(Q,go))/2^dpg n d P go] 
xe X p{i[S n (Q)-S (Q)]/X}. 



(10) 



We then change from integration along C to integration 
along the semicircle C arc at \Q\ — > oo, Im Q > 0, and 
contour C that for /i < goes above the real axis from 
— oo to Q — +0 around the branch cut on the imaginary 
axis, see Fig. [2] (a). We use that the classical trajecto- 
ries Q(t; g) for the Hamiltonian g(Q, P) are expressed in 
terms of the Jacobi elliptic functions [5]. For each g, this 
expression provides conformal mapping of the half-plane 
Im Q > (with a branch cut) onto a g-dependent region 
on the plane of complex time r. We define t(Q, g) as the 
duration of classical motion from the turning point ap(g) 
to Q. Then, for /i < the region on the r-plane that cor- 
responds to the half-plane Im Q > is the interior of a 
rectangle shown in Fig. [2] (b). 

Using that r(Q,g) = dS/dg, we write the exponent in 
Eq. pOj as 



L[S n (Q)-S (Q)} = ^ 



dgr{Q,g). 



(11) 



As seen in Fig. [2] (b) , for any Q on contour C arc and any 
Q' on contour C , Im r(Q,g) < Im r(Q',g). Therefore 
the integral along C can be disregarded. 
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FIG. 2. (a) The contour of integration C for calculating the 
matrix element J7J in the WKB approximation and the aux- 
iliary integration contours for fj, < 0; as = as(g n ) is the 
branching point of P(Q,g„), a B (g) = i{g - ffmin - ^) 1/2 , see 
Eq. (pi. (b) Mapping of the half- plane Im Q > (with a 
branch cut) on the interior of a rectangle on the r-plane for 
/i < by function Q(r;g) that describes the classical Hamil- 
tonian trajectory with given g; Tp , 2ro, and r* are the real 
and imaginary periods and the pole of Q(r;g), respectively. 
The solid (C), dashed (C arc ), and dash-dotted (C') lines are 
the maps of the corresponding contours in (a). The arc in the 
upper right corner is the map of the imaginary axis of Q from 
a,B(g) to ioo; tb is the time for reaching as (gn)- 



On contour C arc t(Q, g) is given by the position of the 
pole t*(<?) of function Q(r;g) [8 , whereas from Eqs. ^ 
and (B|) and the expression for hi(Q,—P(Q,go)) (see 
Supplemental Material) the prefactor in h + is oc 1/Q. 
Then from Eq. 



(ipnlhilipo) ~ o c « c o ex P 



-A" 



dg Im r* (g) 



.'/i i 



(12) 



'/u 



Equation ( 12 ) gives the matrix elements of the fast- 
oscillating field h\ in the explicit form, including both 
the exponent and the prefactor. An excellent agreement 
of Eq. ( 12 ) with numerical calculations is seen in Fig. [3] 
The matrix elements of h 2 are exponentially smaller than 
those of hi and can be disregarded. 

Expression ( 12 ) determines the matrix element h ICS = 
| (ip n \hi |V'o)/v / 2 in Eq. (|6|) for the interwell oscillation fre- 
quency One should compare the exponent in Eq. ( 12 ) 
with the RWA tunneling exponent log 6 go. The latter 
is determined by action So(—a,R(go)) for moving un- 
der the barrier from one well of g(Q, P) to the other 
8J. By symmetry, it is given by twice the real part of 
Eq. (fn) for g n = and Q = +0. From Fig. |] (b), 
Im T\Q,g) — Im To(g)/2 > Im T*(g) for g < 0. Therefore 
for g n < not only |(^„|/n|^ )|, but also \(i/j n \hi\tfj )\ 2 
are exponentially larger than 5 go- As shown in Supple- 
mental Material, this relation holds also for /i > 0. 



IV. DISCUSSION OF RESULTS 

Of utmost interest for observing non-RWA tunneling is 
resonance of hup with states close to the barrier top in 
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FIG. 3. Left panel: the scaled tunnel splitting i>_ = 
(ei — eo)/A for the quasienergy states that maximally over- 
lap with the lowest-;? states (i/'o =t W )/\^ f° r A* ~ 0-95 
and A = 0.08. The quasienergies eo,i are obtained numeri- 
cally from the Schrodinger equation with the full Hamiltonian 
H(t). The dotted lines show a comparison of the peak shapes 
with Eq. (12 1 for renormalized g„, h rcs ; the resonating g n are 
near th e b arrier top in Fig. [T] Right panel: A comparison 



of Eq. (12 1 for h„o — {ipn\hi\tj}o) calculated as a continuous 
function of g (the solid lines) with numerical calculations; the 
triangles and crosses refer to symmetric and anti-symmetric 
states tf> n - A narrow vicinity of g n — where expression ( 12 1 
logarithmically goes to zero is not shown. 



Fig.[l] This is because the matrix elements ( 12 ) fall down 
exponentially with increasing g n , whereas the condition 
of resonance with only one state n requires that \g n — 
g n ±i\ ^> hies', this condition is violated deep inside the 
wells of g(Q, P), because tunnel splitting quickly becomes 
smaller than h rcs with decreasing g n . 

For hjjp ~ At/ one should take into account tran- 
sitions Vo T ~~* V'n v i a intermediate states, which ap- 
pear in higher-order in h. From Eq. (12), if the inter- 



mediate states are arranged in the order of increasing 
quasienergies, the resulting transition matrix elements 
have the same exponent as h rcs . Thus the corresponding 
virtual processes just renormalize the prefactor in the 
resonant tunnel splitting compared to Eq. (12). Figure [3] 
shows extremely sharp resonant peaks of the splitting of 
quasienergy levels for hujp resonating with the RWA in- 
ter/level distance (F 2 /6j)(g n —go)', this distance is also 
renormalized in higher orders in h and was used as an 
adjustable parameter. 

Even moderately weak relaxation modifies the inter- 
well transitions if the oscillator decay rate T exceeds the 
tunneling frequency w s i^_. We will consider the resonant 
case, where the dimensionless decay rate k = T/ui s \ <C 1 
but k ^$> A. If k >• /ires/A, resonant transitions ip Q ' r ' — ¥ 
ijj n occur at dimensionless rate ~ h 2 es /X 2 K. From the 
resonantly excited state ip n the system drifts down in 
quasienergy and approaches the states ^q and i/jq' with 
equal probabilities. As a result, instead of tunneling the 
system incoherently switches between the wells with rate 
~ h 2 cs /X 2 K, see Supplemental Material. 

Relaxation leads to intcrwell switching on its own 
via the mechanism of quantum activation [8] 110] . The 
dimensionless switching rate is vqa ~ texp(— Ra/X). 
From the explicit form of the activation exponent Ra 
[S] and Eq. Q it follows that, for T = 0, the rate of 



/i-field induced switching is exponentially higher than 
the quantum activation rate for fi > —0.35, if the res- 
onant quasienergy level is near the barrier top, \g n \ < A. 
For T exceeding a small /i- and K-dcpcndcnt threshold 
value (still T -C hojF/2ks), Ra becomes smaller than 
the leading-order term in 2A| log h res \. The difference be- 
tween these quantities quickly falls down with increasing 
//, and for realistically not too small A and small quan- 
tum activation prefactor k, the ft,-field induced switching 
may still dominate at resonance. 

In conclusion, we have found a new mechanism of tran- 
sitions between coexisting vibrational states of a para- 
metric oscillator. The transitions correspond to reso- 
nant tunneling in a A-type configuration of quasienergy 
states and come from the terms, which are convention- 
ally disregarded in the RWA. Using the conformal map- 
ping technique, we show that the transition amplitude 
is exponentially larger than the RWA tunneling ampli- 
tude. It displays sharp resonant peaks as a function of 
the modulation frequency. These peaks should make it 
possible to observe the effect in experiments on oscilla- 
tors with a high quality factor, including the currently 
studied Josephson junction based oscillators [I]. 

The research of VP and MID was supported in part by 
the NSF, grant EMT/QIS 082985 



Appendix A: Classical trajectories and conformal 
mapping 

The classical solutions Q(r,g) for g < in terms of 
the Jacobi elliptic functions are given in Ref. Sj except 
that here we count time off from the larger- rather than 
the smaller-Q turning point. In the main text we ana- 
lyzed the conformal mapping of the half-plane Im Q > 
and the integration contour in Fig. 2 (a) onto the r- 
plane for the case /i < 0. The mapping changes for 
g < —(1 — fi) 2 /4, in which case the momentum branching 
point asig) = i(g — ffmin — m) is on the real Q-axis. 
For the ground state, g = go w g mm , this happens for 
fi > 0. On the other hand, for the resonant quasienergy 
g n m g min + hhjp/E s \ close to the barrier top, as (<?«.) is 
imaginary for all |/i| < 1. The mapping depends also 
on the sign of the parameter mj of the Jacobi elliptic 
functions, which itself depends on /i and g. 

The mapping for g < and /i > can be found in the 
same way as for /i < 0. To find the mapping for g > 
we use the expression for the Hamiltonian trajectory 
Q(t; g), which can be obtained by analytically continuing 
the result for g < 0, 



Q(r;g) = 



2 3 /V /2 dn(2(l+i) 5 1 /4 r ) 
k + - K-cn (2(1 + Vjg 1 /^) 



where k± = (l + /i ± 2ig 



1/2 



1/2 



and 



to j (. 9 ) = 1/2 + z(1-/i 2 -4 5 )/8 v ^. 



(Al) 



(A2) 




-7-W/2 o| R^ 
g < 0, m,j < 




on C is exponentially smaller than the integrand on C a 



5 < 0, mj > 




-7-W/4 0| Re? 

< ,g < 1 + ff min 




3 > 1 + Smii, 



FIG. 4. Conformal mapping of the half-plane Im Q > (with 
the branch cut) on the r-plane using the classical Hamiltonian 
trajectories Q(t; g) for underbarrier (g < 0) and overbarrier 
(g > 0) trajectories; in (a) and (b) /i < and /i > 0, re- 
spectively; in (c) < g < 1 + g min , and in (d) g > 1 + p min . 

Parameters Tp and 2to are the real and complex periods of 
Q(r\g), and t* is the smallest in the absolute value pole of 
Q(r;g). The solid line C, the dashed line Care, and the dash- 
dotted line C' show the mapping of the corresponding integra- 
tion contours, see Fig. 2 (a) of the main text. The green line 
is the map of the imaginary Q-axis from the branching point 
a s(fl) to ioo; tb is the time for reaching bs(j„). The blue 
line corresponds to moving from Q — — to Q = +0 along 
the imaginary Q-axis around the branching point as (g) in (a) 
and around a,B(g) and the turning point P(Q;g) — in (c) 
and (d) 



yu 



dgImr(Q axc ,g)< rfjImr(Q',j) (A4) 

J go 



for any Qarc on contour C arc and any Q' on contour C . 

From inspection of Fig. [I] (a), (b) we see that the 
inequality holds for g n < 0, since Im T(Q alc ,g) < 
Im t(Q', g) for all Q arc and Q' . The latter inequality fails 
for g > in a certain range of Q' . However, we emphasize 
that Eq. ( A4 ) is only a sufficient condition for disregard- 
ing the integral along C, and condition Im T(Q aj:cl g) < 
Im r(Q',g) is only a sufficient condition for Eq. (|A4[) to 



hold. Clearly, the matrix element (V'nl^ilV'o) is given by 
the integral along C arc for small g n > 0, where the major 
contribution to the integrals (A4) comes from the range 
g < 0. The range of small \g n \ is of utmost interest for 
observing the non-RWA switching, as explained in the 
main text. 

We note that in Eq. (10) of the main text we dis- 
regarded the change of ^q{Q) for Q > due to the 
tail of the part of this wave function localized in the 
Q < 0-well of g(Q,P). In other words, we approxi- 

mated MQ) = [#(0) ± 4\Q)}/^ by $\Q)/>/2 
for Q > . This is certainly justified in the case of inter- 
est where the matrix element (tp n \hi\ipo) is exponentially 
large compared to the matrix element of tunneling be- 
tween the wells of g(Q, P) for g = go- 

Our numerical analysis shows that the major contri- 
bution to the matrix element (VVil^ilV'o) comes from in- 
tegration along C arc even for large g n > 0. This is is a 
consequence of fast oscillations of the integrand on the 
contour C 



Appendix B: The prefactor in the WKB matrix 
elements 



From Eq. (All, the smallest in the absolute value position 



of the pole of Q(r,g) in the upper half-plane of r is 



,{g)=i 2 



-3/2 



-^cn- 1 



(«+/«-) 



(A3) 



Figure [4] shows the mapping onto the r-plane of the 
integration contours C, C and C arc on the Q-plane plotted 
in Fig. 2 of the main text. It follows directly from the 
explicit expressions for Q(r;g). Panels (a)-(d) refer to 
the four different possible situations. In all cases except 
for the one in Fig. H (b) , the half-plane Im Q > , with 
the branch cut along the imaginary axis, is mapped onto 
the interior of a rectangle. In Fig. [4] (b), the upper part 
of the rectangle is replaced with a curve that corresponds 
to going between the branching points of P(Q, g) on the 
real Q-axis. Panel (a) coincides with Fig. 2(b) of the 
main text and is given for completeness. 

The integral along C can be disregarded when calcu- 
lating the matrix elements (V'nl^ilV'o), if the integrand 



It follows from Eqs. (10) and (11) of the main text and 
the ensuing analysis that, in the WKB approximation, 
the matrix elements of operators hi t 2 are given by inte- 
grals over an arc on the Q-plane that goes from to 7r 

for |Q| — > oo, 



{ip n \hj\ipo} « c n c Re / dQhj(Q,-P(Q,g )) 
Jc arc 
r i rsn 
x (dpg n d P g )~ 1/2 exp 



A 



dgr{Q,g) 



.'/(i 



(Bl) 



where j = 1,2. The explicit form of hi t 2 is 

k = (P - iQ) [{Q 2 + P 2 )/6 - m/4] (P - iQ) 

+ (Q 2 + P 2 )/2, 
h» = (P~ «<2) 4 /24 + (P - iQ) 2 /4. (B2) 



The normalization constant c n in Eq. (Bl| is c ri 
[ri 1) ( 5n )/4]-V2 i for Sn >0. 



On the integration contour in Eq. (Bl ), from Eqs. (3) 
and (9) of the main text we have P(Q,g) ps —iQ + L 
d P g ps -2iQ 2 , T(Q 7 g) ps T *(g). Then from Eq. flB2J ) 
fti(Q,-P(Q,5o)) « 2Q/3 and h 2 (Q,-P(Q,g )) ps 
—5/24. Substituting these expressions into Eq. ( |B1[ ) gives 
Eq. (12) of the main text for the mat rix element of h±. 
In the case of h 2 the integrand in Eq. (Bl I is ex 1/Q 2 for 



\Q\ — > oo, so that the integral (Bl) goes to zero. There- 
fore the matrix element (tp n \h2\1po) is determined by the 
exponentially smaller contribution that comes from inte- 
gration along the path C. 



Appendix C: Effects of dissipation 

In the presence of dissipation, the dynamics of the os- 
cillator in the rotating frame is described by the master 
equation for the density matrix p, 



d T p = iX 



P,9 



up, 



kp = n(n + 1) 



x (a* ap — 2apa* + pa' a) + n(aa' p — 2a' pa + paa '()C1) 

Here, k ~ T/uj s \ is the scaled coefficient of viscous fric- 
tion (the classical friction force is —2Tdq/dt) and n — 
[exp(hojp/2kBT) — 1] is the oscillator Planck num- 
ber; a) and a are the raising and lowering operators, 

a = (2X)- 1 t 2 {Q + iP). 

We will consider /i-induced switching between the 
period-two vibrational states in the most interesting case 
of moderately small decay rate, 1 > k > h lcs . For 
resonant transitions "00 ~^ ^n the detuning A — 
X~ 1 (g n — go) — (luf/lo s \) will be assumed small, |A| <C k. 
If the system is initially prepared in one of the wells of 
g(<5,P), over dimensionless time ~ kT 1 the distribution 
p becomes quasi-stationary, except for the well popula- 
tions, which as we show vary on a much longer time. 

For k <C 1 off-diagonal matrix elements of p on 
the eigenf unctions of operator g are small. They can 
be disregarded except for the matrix elements p^ ' — 

(4'n\p\4'o ) f° r the resonant state ip n and their conju- 
gate (the other exception are matrix elements in a narrow 
range of <?-values near the barrier top where the tunnel 
splitting is ~ Xk). From Eq. (CI I, to the leading order 
in /i res /K, in the quasi-stationary regime 



~(r) 

PnO 



h ICS /X 



k(W + W n ) 



( ~(rr) 

[Poo 



(C2) 



,-,(') 



and similarly for p n Q. Here we use the interaction rep- 
resentation, p = V'pV with V = V(t) = cxp(—igr/X); 
p nn and /5q are the diagonal matrix elements of p on 

(r) 

the wave functions ip n and "0g , respectively. 

Parameters Wq and W n are determined by the rates of 

(r) 

dissipation-induced transitions from states ip and ip n . 



FromEq. ( |Cl] ), 

w m =y~] W mm <, 

W mm , = (n + l)|(VwH0m)| 2 +nK^Ja|0 ro <}| 2 (C3) 

Matrix elements W mm i were calculated before [5]; they 
exponentially decay with increasing \m — m'\, as is also 
clear from Eq. (12) of the main text. In Eq. (C2) we 



neglected direct tunneling between the states ip$ and 

The matrix elements ( C2 ) determine influx into the 
population of the resonant state p nn . If initially the os- 
cillator is in the right well of g(Q,P), for example, then 

(rr) (ID 

Poo -^ Poo ' 1 Pnn for not too large r, and from Eq.(Cl) 
the influx term is 



l h Tm n (r) ss 2vu?> {rr) 

v h = X- 2 h 2 J{ K (W + W n )}. 



[d T pnn]m — —2X 



(C4) 



Dissipation leads to transitions from the resonant state 
ip n to neighboring states ip m , with rates 2KW nm . They 
result in a current in the space of state populations. The 
current density can be found from the balance equa- 
tion for the diagonal matrix elements p mm assuming 
that Pmm is a smooth function of to, i.e., p m +k,m+k ~ 

Pmm + kdmPmrn for |fc| < TO. 

In the WKB approximation, the matrix elements 
Wmm' for small |to — m'\ are expressed in terms of the 
Fourier components of the periodically oscillating in time 
momentum P and coordinate Q on a classical trajectory 
with Hamiltonian g(Q,P) [T5] . It is straightforward to 
show that, for a smooth p mm . 



II! ' 



(W m ' m p. 



m. m f J m' rrr 



Wmm' Pmm) « A 1 9 m [I(m)p n 



(C5) 

where I(m) is the classical action for an orbit with 
g(Q,P) — g m ', from the Bohr quantization condition, 
d m I(m) = X. 

Equation d T p m m — — (2K/X)d m [I{m)p mm ] has a stan- 
dard form of a continuity equation. It applies for g m 
not too close to the resonant quasienergy g n , the barrier 
top, and the minima of g(Q,P). Inside the wells, for 
\g m \ S> X and for the tunneling splitting small compared 
to Xk, subscript to enumerates intrawell states. For such 
states, in the time range kT 1 <C t <C n(X/h res ) 2 the prob- 
ability current I(m)p mm is constant. Since the wells of 
g are symmetric and there is no accumulation away from 
the minima of g, the current is equal to half the influx 
into the resonant state (C4). Since the current into the 



initially empty well gives the rate of interwcll switching, 
this rate is equal to v^. 

We note that, even for zero temperature T, dissipa- 
tion leads to population of excited eigenstates of g and 
to interwell transitions via quantum activation. For 
T = the quantum activation exponent is Ra — 
21m f dg[To(jg) — 2r*(g 1 )] [8]. Comparing this integral 
with the exponent of v% oc h 2 as we obtained the range of 



/j, where switching via a resonant state is exponentially more likely than quantum activation. 
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